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ABSTRACT 

We measure the ellipticity of isolated clusters of galaxies in the Sloan Digital Sky Survey (SDSS) 
using gravitational lensing. We stack the clusters, rotating so that the major axes of the ellipses 
determined by the positions of cluster member galaxies are aligned. We exclude the signal from the 
central 0.5 h~ l Mpc to avoid problems with stacking alignment and cluster member contamination. 
We fit an elliptical NFW profile and find a projected, two-dimensional axis ratio for the dark matter 
of / = b/a = 0.48^ 09 an d rule out / = 1 at 99.6 per cent confidence thus ruling out a spherical 

halo. We find that the ellipticity of the cluster galaxy distribution is consistent with being equal to 
the dark matter ellipticity. The results are similar if we change the isolation criterion by 50 per cent 
in either direction. 

Subject headings: cosmology: dark matter — cosmology: galaxy clusters — cosmology: large-scale 
structure of universe — cosmology: observations — galaxies: clusters: general 
galaxies: halos — 
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1. INTRODUCTION 

Cosmological simulations can be used to predict many 
different statistics of the mass distribution in the Uni- 
verse. The most commonly employed statistic is the two- 
point correlation function, or its Fourier counterpart the 
power spectrum. Three-point statistics are much harder 
to predict and measure, and higher orders are rarely 
discussed. A more popular statistic is the number of 
peaks in the mass distribution, as given by the number 
of clusters of galaxies. The dark matter power spectrum 
and the number of clusters of galaxies are often cited as 
among the best ways to constrai n the properties of dark 
energy (e.g. lAlbrecht et aLll2006l ). 

Uncertainties on cosmological parameters are de- 
creased when two measurements have different param- 
eter degeneracies and are often referred to as 'comple- 
mentary'. In this paper we consider a statistic which 
may offer complementary constraints on cosmology: the 
shapes of peaks in the mass distribution, as probed by 
the ellipticity of galaxy cluster dark matter halos. In 
addition this may place important constraints on modi- 
fications to the law of gravity since we may compare the 
results from both dark and light matter, as also tested 
by studying the dark and light mat ter distributions in 
the bullet cluster (|Clowe et al.ll2006f ). 

Predictions of clust er ellipticities come mostly from 
numerical simulations ifWest et a l. 1989; dc Thciie ~t al.l 
1995 1 : iJing fc Sutj l2002t jFloor et all 120031: iHo fc White! 
2004 iFlores et all 120051 : iRahman et all l2006f ) . The 
ellipticity is expected t o dep e nd on cosmo l ogica l 
parameters (lEvrard et all 119931: iSplinter et al.l [19971 : 
iBuote fc Xulll997l : ISuwa et al.ll2003l:lRahman et al ll2004f) 
and to evolve wit h redshift (|Kasun fc EvrarcF 120051 : 
lAllgood et al.ll2006f ) an evolution which itself might de - 
pend on cosmology ([Hopkins et alJ[2005t IHo et al.ll2006f ). 
The distribution of sub-halos within a cluster halo is 
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found to be an i ndicator of the overall halo ellipticity 
(|Bode et al.ll2007t ) but is slightly less elliptical. 

The ellipticity of the brightest cluster galaxy and the 
ellipticity of the distribution of cluster member galax- 
ies is much easier to observe than any other cluster 
ellipticity mea sure. This has been studied by a num- 
ber of authors (IWest fc Bothun|[T990HPlionis et alj |1991t 



Rhee et al.lll991t IStrazzullo et al.ll2005l ). Most recently. 
Wang et al. studied groups of galaxies in SDSS 



and found an alignment between the brightest cluster 
galaxy (BCG) and the distribution of galaxies which was 
strongest in the most massive groups and between red 
BCGs and red group member galaxies. 

In this paper we focus on measuring the dark matter 
ellipticity directly using gravitational lensing. We also 
compare this ellipticity to the ellipticity of the cluster 
member galaxy distribution, to sec how reliable a tracer 
of ellipticity the light is, and therefore whether it can be 
used by itself for cosmological studies. 

Gravitational lensing has been used very successfully 
to measure the mass and profile of clusters of galaxies 
by many authors. The work mo s t relevant to o u r study 
is t hat of ISheldon et all rf200l : ISheldon et all (l2007bJ ) 
and lSheldon et al.l ( 2007al) who stack the lensing signal 



from many clusters of galaxie s to find an average signal. 
iNatarajan fc R efrcgicr (2000) proposed that when stack- 
ing the shear signal from many halos the stack could be 
made while retaining information about the major axis 
of the halo, as observed from the distribution of light. 
The stacked shear map should then provide a constraint 
on the ellipticity of the dark matter halo, if indeed the 
mass and light were aligned. They considered lensing by 
galaxies but we apply the same technique to clusters of 
galaxies here. 

The ellipticity of the dark matter distribution has 
been studied observationally for galaxy sized halos using 
gravitational le nsing using data from the Red-Sequence 
Cluster Surv ey feoekstra et "all 1 2004) . the Sloan Digital 
Sky Survey (jMandelbaum et all 12000 ) and the Canada- 
Franc e-Hawaii Telescope Legacy Survey (jParker et al.1 
l2007f h This has proved to be extremely difficult and 



2 



we present here, for t he first time, results fro m stacking 
galaxy cluster halos. ICvpriano et al.1 (|2004f ) has previ- 
ously made measurements of the cluster dark matter el- 
lipticity using gravitational lensing measurements from 
individual clusters. They found a good agreement be- 
tween the dark matter halo orientation and the orienta- 
tion of the brightest cluster galaxy. 

When calculating angular diameter distances and the 
mean density of the Universe we assume a flat cosmology 
with Q m = 0.3. The NFW halo profile has a very weak 
dependence on the fluctuation amplitude. We assume 
erg = 0.8 for this calculation. 

The structure of this paper is as follows. In Sec. we 
describe our data set and operations we have performed 
on the data such as rotating and stacking and removal 
of neighboring clusters. In Sec. [3] we use two theoreti- 
cal models - a Singular Isothermal Ellipsoid (SIE) and 
a Navarro-Frenk- White (NFW, see iNavarro et 111111997?) 
model - of the mass and light distribution, and look at 
correction factors from the redshift distribution and clus- 
ter decontamination. Our results, and some interpreta- 
tions of these, are presented in Sec. |U Conclusions are 
summarized and discussed in Sec. 

2. DATA 

In this section we describe the catalogues used and the 
operations we carried out before comparing with models, 
including cluster selection and stacking and rotating. 

2.1. Catalogues 

We use the cluster catalogue of Ko ester et al.l (|2007aD . 
which is at the present time the largest existing galaxy 
cluster catalogue, consisting of 13,8 23 galaxy cluster s 
from the Sloan Digital Sky Survey (|York et al.l I2000D . 
The cluster galaxies are red-sequence members (occu- 
pying the so-called E/SO ridgeline in colour-magnitude 
space), brighter than 0.4L* in the i band and between 
rcdshifts 0.1 < z < 0.3. They also lie within a circular 
aperture of radius i?200, which is the estimated radius 
within which the density of galaxies with 24 < M r < 16 is 
200 times the mean density of such galaxies. The number 
of cluster members inside this aperture is 10 < < 
188; the lower limit is a r equirement for inclus ion into 
the catalogue. We refer to iKoester et al.l (|2007al) f or de- 
tails of the catalogue and to IKoester et al.l (|2007bf ) for a 
description of the cluster selection algorithm. 

In order to define an isolated sample, we removed clus- 
ters found to be too close to each other, as seen on the 
sky. Details on this close neighbour removal can be found 
in Sec. 12.21 Clusters that are too close to the survey 
edge are also removed, with the requirement that the 
minimum distance from the cluster center to the survey 
edge is Ihr 1 Mpc. We are left with a total of 4281 clus- 
ters to analyse for our purposes. We have also organised 
the clusters into 4 redshift bins between z — 0.10 and 
z = 0.30, each of width 0.05 in z. Since these cluster 
rcdshifts are photometric, with uncertainties of ~ 0.01 
(Koester et al 2007a), the spectroscopic bin width will 
be very slightly larger than the photometric bin width. 
We do not take into account this small broadening in our 
analysis because we expect it to have negligible effect on 
the ellipticity results. 

The shear galaxy catalogue is the same as used in 
ISheldon et al.l (|2004l ). except that the area covered is 



larger, ~ 6325 deg . There is approximately one galaxy 
per square arcminute in this catalogue. Galaxies in the 
shear catalogue have extinction-corrected r-band Pet- 
rosian magnitudes less than 22. Stars have been re- 
moved fro m the catalogue by th e Bayesian method dis- 
cussed in iScranton et alj ([2002) . Unresolved galaxies 
and objects with photo-z errors greater than 0.4 have 
also been removed. To correct the shapes of galaxies for 
eff ects of PSF dilution and anisotropy, the techniques 
of iBernstein fe Jarvisl 1 2002) were used, with modifica- 
tions specified in Hirata fe SelialJ ((2003). We refer to 
ISheldon et al.l (|2004l ) for full details of the compilation 
of the shear catalogue. 

We do not want to include in the shear catalogue galax- 
ies which are already in the cluster member catalogue. 
Therefore we remove from the shear catalogue all galax- 
ies which are close, as seen on the sky, to a cluster mem- 
ber. We cut on a physical (as opposed to angular) dis- 
tance of 0.012 hr 1 Mpc, as calculated at the redshift of 
the cluster. This distance corresponds to 5 arcseconds 
at redshift z — 0.2. In principle this cut may also re- 
move some real background galaxies, however the shear 
measurements of these galaxies will in any case likely be 
adversely affected by the light contamination by the clus- 
ter members. Note that we do not eliminate all cluster 
members from the shear catalogue with this cut, only 
those clusters that fit the selection criteria set by the 
maxB CG selection method (red galaxies, brighter than 
0.4 L*; IKoester et alj[2007ah . 



2.2. Postage stamp size and close neighbor removal 

In order to avoid contamination of the shear signal 
from neighboring structures, we only include in our anal- 
ysis galaxies which are sufficiently close to a cluster, as 
seen on the sky. We decide this distance by consider- 
ing the predicted contribution of neighboring clusters to 
the shear signal, and then use this to decide the required 
separation of neighboring clusters. 

We include in our shear catalogue galaxies in a square 
postage stamp of 10 x 10 h^ 1 Mpc center ed on each clus- 
ter. T his choice was made using Fig. 8 of lJohnston et al.1 
(|2007f ). which shows model fits for the lensing signal 
split up into contributions from several components. We 
want the lensing signal in our postage stamp to be domi- 
nated by the central cluster, not by the contribution from 
neighboring mass concentrations such as nearby clusters 
and filaments. Our interest therefore lies in comparing 
their NFW profile (green line) to the contribution from 
neighboring halos (blue line). The results in their Fig. 8 
are shown for 12 richness (N^®) bins. The mean of 
the Ngl° B ° values of our cluster sample is ~ 24, directing 
us to the panel N^^[18 — 25]. We match our postage 
stamp limit at the radius where the NFW contribution 
is roughly equal to that of the contribution from neigh- 
boring clusters. Note that we are also removing close 
neighbors (see below), so the contribution from neigh- 
boring halos would actually be even lower for our case, 
so our cut may be conservative. 

Clusters that are too close to each other in projected 
separation represent a challenge for our analysis. The 
mass distribution of a cluster will affect the shear field 
of its neighbor. This could have a significant effect on 
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the measured ellipticity. While it is possible to isolate 
clusters in three dimensions within a simulation, we are 
plagued by overlaps on the sky since gravitational lcnsing 
measures the projected mass. We therefore endeavor to 
remove this complication by selecting relatively isolated 
clusters for this analysis. For proper comparison with 
theory, a similar sample should be made from simula- 
tions. However, in this paper we are primarily concerned 
with the first significant detection of ellipticity in a large 
sample. 

We do not want the shear from a neighbor to appear 
within our postage stamp, but this is unavoidable to 
some extent, because the shear at one point is affected 
by mass far away (shear is non-local). The effect on the 
measured ellipticity of having neighbors is twofold: (i) It 
can make the distribution more circular, if the neighbor 
position is uncorrelated with cluster major axis. This 
will occur due to chance alignments close to the line 
of sight, particularly from clusters in different redshift 
bins; (ii) It can make the distribution more elliptical, if 
the neighbor position is correlated with major axis di- 
rection. Pairs of clusters which are physically close will 
prob ably be aligned alo ng the major axis of the cluster 
(e.g. iPlionis et al.lll991h . 

We therefore remove clusters that are closer than bh^ 1 
Mpc (corresponding to half the width of the postage 
stamp) to neighboring clusters. When a cluster has 
one or more neighbors with an angular separation corre- 
sponding to less than 5/i _1 Mpc, calculated at the middle 
of the redshift bin of the cluster (see later), we discard it 
if any of the neighbors have a higher N^®® . The neigh- 
bor removal was done consecutively from low to high 
redshift. This process reduced the number of available 
clusters to 4281. In a later section, we analyze the use of 
different minimum close neighbor distances. 

Although we remove a cluster's less rich neighbors from 
our sample, the shear pattern of the remaining cluster 
will already have been affected by its neighbor (s). How- 
ever, by retaining the richest of the neighboring clusters 
we hope that the shear field is dominated by this cluster. 

2.3. Stacking and rotating 

On average, we have only around 1 shear galaxy per 
square arcminute and the uncertainty on the shear for 
a single galaxy is an order of magnitude larger than the 
shear we are trying to measure. Therefore we need to use 
the shear signal from many clusters in order to obtain a 
significant signal. We therefore stack the clusters on top 
of each other to improve the signal-to-noise ratio. In 
other words, we use information from the postage stamp 
field of shear galaxies for all clusters simultaneously. 

The stacking could be carried out in either physical or 
angular space. For our method the two approaches are 
exactly equivalent if the redshift bins are small enough. 
We stack in angular space and use redshift bins of red- 
shift width 0.05. This causes a radial blurring because 
two clusters of the same physical size will be stacked on 
top of each other in angular space to have different an- 
gular sizes. The blurring of the shear and light maps will 
be of at worst plus and minus 20 per cent (for the low- 
est redshift bin). Since superposing ellipses of different 
scalings retains the original ellipticity, this results in a 
slightly smoother cluster profile but will not affect our 
ellipticity results. We find that our results are fairly sim- 



ilar even when comparing two very different profiles, see 
Sec. El 

Straightforward stacking of elliptical clusters with ran- 
dom orientations would erase any ellipticity and produce 
a circular average cluster. Before stacking, we there- 
fore rotate each cluster to lie along an x-axis, which is 
the major axis as defined by the ellipticity of the clus- 
ter members, s ee Fig. [1] This rotate-and-sta ck method 
is described bv lNatarajan fc Refregierl (|2006) for use in 
galaxy-galaxy lensing, and we have, for the first time, 
applied this technique for use on cluster lensing. 

We calculate the direction of a cluster major axis from 
the positions of the cluster members, as defined in the 
cluster catalogue. We do not take into account the lu- 
minosity of each cluster member. The cluster center 
(%o He) was taken to be the position of the brightest 
cluster galaxy (BCG) as defi ned by the maxBCG algo- 
rithm (|Koester et al.|[2007H ). However, the position of 
the BCG is not necessarily coincident with the cluster's 
actual centre of mass. For comparison, we therefore cal- 
culate the center of each cluster as given by the mean 
position of the cluster members. The mean physical off- 
set (for redshift bin 0.20 < z < 0.25) between the two 
center definitions is ~ O.lbh^ 1 Mpc, and the standard 
deviation 0.09/i _1 Mpc. Compared to our mask radius 
of 0.5/i _1 Mpc, therefore, the shift in centre position is 
relatively small. Any effect this mis-centering does have 
will increase the ellipticity of the members, which we do 
not focus on here, and cause the misalignment angle to 
tend towards the direction from the BCG to the center 
of the cluster member distribution. This would in itself 
be an alternative and potentially useful way to stack the 
clusters to obtain the results presented here. Therefore 
we do not consider this effect further. 

To find the ellipticity angle of rotation of the cluster, 
we use the quadrupolc moments of the cluster members. 
The quadrupole moments are given by: 

Qxx = ({xi -x c ) 2 )i (1) 
Qxy = ((xi - x c ) (yi - y c ))i (2) 

Qyy = {{Vi ~ Vof)i (3) 
where the summation i is over the cluster members. We 
convert this into the ellipticity components e\ and ei of 
the cluster through the relations: 




The angle of the cluster, anticlockwise from positive x 
axis, is then: 

<T* = \ atan(|) . (6) 

We rotate the positions of the cluster member galaxies 
and the positions of the shear galaxies using the following 
transformation 

x Iot =d cos(6> - rot ) (7) 

y rot = dsin(0-0 rot ) (8) 

where (d, 9) are the polar coordinates of the galaxy to be 
rotated, relative to the cluster center. For clusters with 
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Fig. 1. — Total number of galaxies per unit area for redshift bin 3 (0.20 < z < 0.25) for illustration. From left to right: shear galaxies, 
cluster members and total galaxies (shear galaxies plus cluster members). The left and right panels include galaxies which are uncorrelated 
with the cluster. We draw our contours relative to the constant background level (which is the same number in both left and right panels, 
since the cluster member catalogue does not contribute). For the left and right panels, contours are equally spaced from 1.2 to 2.6 times 
the background level, in steps of 0.2. The dashed contour shows 2 times the background level. For the central panel the contours are 0.2, 
0.4, 0.6, 0.8 times the same background level. The dashed circle shows location of mask, see text. For comparison, a 3 X 10 14 h~ 1 Mq 
cluster at a redshift of 0.2 has R200 ~ lAh~ 1 Mpc, where -R200 is the radius within which the mean density of the cluster is 200 times the 
mean matter density at the cluster redshift. We have zoomed in to show the central part of the postage stamp for clarity. 

member catalogue galaxies removed, the middle panel 
Q2 shows cluster galaxies from the cluster member catalogue 

and the right panel shows the sum of shear galaxies and 
cluster members. 

0.01 5 The left hand panel clearly contains a significant num- 

ber of cluster members, despite the fact that the galaxies 
from the clu ster member c a talogue are not included. As 
0.01 ^ described in lKoester et al.l (|2007b| ), cluster members for 

this catalogue were identified using the maxBCG algo- 
rithm. This algorithm employs the red sequence method, 
0-005 based on the observational fact that cluster galaxies oc- 

cupy a narrow region (a so-called ridgeline) in a colour- 
magnitude diagram. This method is designed to conser- 
vatively select red galaxies at the central area of a cluster. 

As an illustration we investigate the central region of the 
stacked cluster in the third redshift bin (0.2 < z < 0.25). 
In the central square arcminute the number of galax- 
ies in the member catalogue is ~ 40% of the number of 
galaxies in the (members removed) shear catalogue, af- 
ter subtracting the constant background level i.e. most 
of the cluster members are not in the cluster member 
catalogue. 

The existence of these extra members allows a very 
convenient check on our stacking and rotating: the ro- 
tation angles were calculated from the cluster member 
catalogue alone, whereas the left hand panel does not 
include the galaxies used to decide the rotation angle. 
Therefore the fact that we see ellipticity in this panel 
means that the angle calculated from the members cat- 
alogue is correlated with the angle of the extra cluster 
members. If, for example, we had calculated rotation 
angles from a very small number of galaxies from the 
cluster member catalogue, there would be a large degree 
of randomness due to shot noise, and the alignment of 
measured and true ellipticity would be random to a large 
extent, resulting in a circular pattern in the left panel of 
Fig.ffl 

The object of this paper is to calculate the ellipticity of 




Fig. 2. — Tangential shear, 7t, as measured from all galaxies 
in the background galaxy catalogue for a zoomcd-in view of the 
postage stamp. To reduce the noise we have smoothed with a 
Gaussian with standard deviation Xh~ 1 Mpc. The location of the 
mask is shown by the dashed circle. 

ellipticities close to zero, 9 Iot in Eq. [6] has little mean- 
ing. From the cluster members alone, 18% of our clus- 
ters have an ellipticity (e = y'ef + ef , using Eqs. [4] and 
[5]) of less than 0.1 and only 5% have ellipticities less than 
0.05. Therefore the angle is reasonably well defined. 2% 
of the clusters have an elli pticity greater than .5. The 
cluster selection criteria bv lKoester et al.l (|2007al) and/or 
our isolation criteria have therefore done a reasonable job 
of identifying isolated clusters. In addition the clusters 
seem relatively undisturbed, i.e. have low ellipticity. This 
also illustrates that it would not be particularly useful to 
bin the clusters according to the cluster member elliptic- 
ity because the range is relatively small (and our final 
signal to noise is quite low). 

Fig. [TJ shows the number of galaxies per unit area, for z 
bin 3, after stacking and rotating. The left panel shows 
galaxies from the shear galaxy catalogue, with cluster 
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the dark matter, as measured from gravitational lensing, 
and compare it with that of the light-emitting galaxies. 
Any misalignment in the stacking and rotating will tend 
to make the dark matter appear less elliptical. However 
this misalignment will have the same effect on the cllip- 
ticity of the light, as measured from the extra cluster 
members (left hand panel) alone. Therefore we can com- 
pare like with like, and assess the relative ellipticity of 
the dark and light matter, despite any misalignment. 

We have chosen the contour levels so that a fair com- 
parison can be made between the left and central panels: 
the outermost contour corresponds to the same cluster 
member density in each. Therefore we see that at large 
radii most of the cluster members are not included in the 
cluster member catalogue. For all the axis ratio mea- 
surements reported in this paper we exclude the central 
regions (see Sec. I3.3|) . this is shown by the dashed circle 
in Fig [TJ Outside this excluded region we see that the 
contours change only a little from the left to the right 
panel. This is convenient because it means that it is not 
too important whether we compare the dark matter el- 
lipticity with the light ellipticity derived from either the 
left or the right hand panel. 

We calculate the tangential, jt, and cross, "fx, com- 
ponents of the shear 7 for each shear galaxy 



7t = 7 cos(2a) 
7x =7 sin(2a) 



(9) 
(10) 



where a is the angle between the shear galaxy major 
axis and a tangent through the center of the shear galaxy, 
with respect to the cluster center. We calculate the shear 
estimate for each galaxy from the polarizations (e) in the 
shear galaxy catalogue 



1 



7 = 



(11) 



where the factor Sgh ~ 0.88 is the a verage respons i yity o f 
the source galaxies to a shear, see iSheldon et ail (|200lh 
and references therein, and e = (a 2 — b 2 ) / (a 2 + b 2 ) , where 
a is the semi-major and b the semi-minor axis of the shear 
galaxy. We have redone our main NFW dark matter 
ellipticity result using a value Ssh which is twenty per 
cent higher, and find only a negligible change. Note that 
the tangential and cross shears are invariant with respect 
to rotation of the cluster coordinates. We bin the shear 
galaxy catalogue into square pixels of size 0.4x0.4 arcmin 
on the sky, but our main results arc insensitive to this 
exact value. We average the shear estimates in each pixel. 

Fig. [2] shows the smoothed tangential shear. As ex- 
pected, the tangential shear is largest in the cluster cen- 
ter. To interpret this figure further it is helpful to con- 
sider the tangential and cross shears pattern predicted 
from popular cluster models, discussed further in Sec- 
tion[3l Cluster ellipticity causes little, if any, cross shear, 
depending on the cluster profile. For an elliptical SIE the 
cross shear is exactly zero. For an elliptical NFW with 
a major to minor axis ratio of 0.5 (our best fit result), 
a mass of 1 x 10 /i _1 M0 and cluster redshift 0.15, the 
maximum cross shear occurs when approaching the cen- 
ter of the lens. Just outside our mask radius of 0.5h~l 
Mpc, the maximum cross shear is 0.0012 for a source red- 
shift of 0.3. This value is 10 per cent of the maximum 
tangential shear outside the mask, and is therefore small 



compared to our uncertainties. The main effect of clus- 
ter ellipticity is to produce an ellipticity in the tangential 
shear map. Therefore the tentative visual indication of 
some horizontal elongation in this figure is our first hint 
of dark matter ellipticity. 

To compare the data with models we must incorporate 
the errors on the shear estimates. The errors on the shear 
measurements are given by 

(12) 



7t 



'SN 



where o~i is the uncertainty in the shape measurement due 
to the finite number of photons falling in each detector 
element, plus detector noise, and ctsn is the 'shape noise' 
due to intrinsic variance in the unlensed galaxy shapes 
(assumed the same for all galaxies). 

We calculate the shape measurement uncertainty <ii 
from the uncertainties in the two components of the el- 
lipticity ei and e 2 using 



1 a e i + CT e2 

25sh 2 



(13) 



The first factor converts from e to 7 and the second part 
assumes that the uncertainty on the two ellipticity com- 
ponents is essentially equal and uncorrelated, which is 
a pproximately true f or some shear estimators (e.g. Fig. 2 
of lBridle~e t al. 2002]). If these assumptions are true then 
it follows that the ellipticity uncertainty on a component 
is independent of rotation. 

We estimate the shape noise asN by calculating the rms 
dispersion in 7r values as a function of <Xj . We compare 
this to (J 7T for various osn values and find ctsn = 0.24 to 
be the best fit. 

To find the error on the shear for each pixel we take 
the mean of the errors ct 7t and divide by the square 
root of the number of galaxies in the pixel. Galaxies will 
have a range of sizes and therefore of measurement er- 
rors. When we calculate the average shear in each pixel 
we do not weight the galaxies according to their elliptic- 
ity errors. This is because this would tend to upweight 
the better measured galaxies. As better measured galax- 
ies may preferentially tend to be cluster members that 
have leaked into the shear galaxy catalogue, this might 
preferentially weight up cluster members. This would 
have to be taken into account when removing the bias 
on the shear due to cluster member contamination (see 
Sec. 13. 3|) . which would be difficult. We therefore use 
all the galaxies in the shear galaxy catalogue without 
weighting. Due to the cuts already made in creating the 
shear galaxy catalogue, the difference in weight between 
the noisiest and least noisy galaxies is only 30% therefore 
the weighting would not make a large difference to our 
analysis. 

3. MODELLING 

3.1. Mass and light distributions 

The shear of galaxies depends on the mass and struc- 
ture of the cluster acting as a lens. To model the clus- 
ter mass distribution, we use two alternative theoreti- 
cal models: a Singular Isothermal Ellipsoid (SIE) and a 
Navarro-Frenk- White (NFW) model. The NFW model 
is preferred from simulations, but we also include results 
from the simpler SIE model to show an extreme and sim- 
ple example of the dependence of our results on the clus- 
ter profile. 
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Fig. 3. — Left panel: Relative probability distribution of the axis ratio /lm = b/a for the galaxy number density map from the shear 
catalogue plus member catalogue, thus including all cluster members. This includes clusters at all redshifts for an NFW model (dark, red 
line) and a SIE model (light, green line). The dashed line shows the result from shear catalogue galaxies only (for an NFW model). Right 
panel: Relative probability distribution for the dark matter axis ratio /dm- The dark, red line shows result from the NFW model, and the 
light, green line shows result from the SIE model. The dotted vertical line shows /dm = 1 corresponding to a circular distribution. 



The SIE model corresponds physically to a distribution 
of self-gravitating particles with a Maxwellian velocity 
distribution with one-dimensional velocity dispersion a v . 
The convergence (normalized mass density) k — E/E cr ; t 
of an SIE is given by 

- 2 D, 



K = 2tt- 



(14) 



where r is the generalized radius 

r=(x 2 / + y 2 r 1 ) 1/2 (15) 

and f — b/a is the axis ratio of the ellipse (b < a), a v 
is the velocity dispersion and c is the speed of light, x 
and y are coordinates in the plane of the sky, at the 
cluster redshift. The distances Z?ds and D s are the angu- 
lar diameter distances between the lens (deflector) and 
the source, and from the observer to the source, respec- 
tively. The SIE peaks sharply in the central parts, but 
we mask out the central regions (see Sec. I3.3|) . For the 
SIE model, we have the simple relation that the normal- 
ized surface density equa l s the tangential shear 7x = n 
(|Kassiola &: Kovnerlll993t iKormann et al.iri994h . 

The NFW model is a more complicated but more real- 
istic model based on numerical simulations. In order to 
implement the NFW, we need the projected mass of the 
cluster. To c alculate the projected mass we use the equa- 
tions g iven in lWright &: Brainerdl 1)20001 ) and I B artelmannl 
(1996) using our generalized radius of Eq. [15] to make 
the cluster mass distribution elliptical. We use M 2 oo, 
the mass enclosed within the radius at which the density 
is 200 times the mean density of the Universe, for con- 
sistency with simulations. We derive the concen tration 
param eter, c, where c oc M@ according to Eq. 12 of lSeliakl 
(2000), where we interpret the virial mass M as -A^oo- 
We use (3 = —0.15, as appropriate for an NFW model. 
The shear for an elliptic al mass distri bution is calculated 
using the eq uations inlKeetonl (|2001h which are derived 
from those in [Schramm (1990). A she ar map using these 
equati ons is illustrated in Fig. 1 of iBridle fc Abdallal 
(2003). 



We calculate probability as a function of our free pa- 
rameters in each redshift bin, and marginalise over all 
but the axis ratio /. We then obtain a single result for 
/ from combining all the redshift bins by multiplying 
the probabilities from different redshift bins together for 
each / value. This is the correct calculation if we believe 
that the other parameters have different values in each 
redshift bin, but that / is the same for all redshift bins. 

3.1.1. Estimation of the light matter axis ratio (J-lm) using 

positions 



When stacking the clusters fSec. 12.3( 1 . we calculated in- 
dividual cluster ellipticities based on the cluster galaxies 
in the members catalogue. We do not use these results as 
our measure for the light matter ellipticity for the stacked 
cluster because they are relatively noisy due to the small 
number of members (~ 10 for the least rich clusters). 
Furthermore, we know that the members catalogue does 
not in fact contain all the cluster members, and has some 
selection criteria that may affect the ellipticity (require- 
ment on proximity to cluster centre). We therefore use 
a x 2 analysis to find the light matter ellipticity of the 
clusters. We model the light map as coming from (i) 
non-cluster galaxies which have a constant density across 
the postage stamp plus (ii) a contribution from the clus- 
ter galaxies which is assumed to have a galaxy density 
proportional to the mass profiles of Eq. [T4l 



i prod (r) = Kk(t) +n Q 



(16) 



where uq is the background level of galaxies per pixel and 
K is a constant. Note that we do not assume that the 
light is some constant multiple of the mass, only that 
the light map is proportional to an SIE or NFW pro- 
file (which may have different parameters than the dark 
matter distribution). We do not tie the dark and light 
map parameters together because we wish to investigate 
whether the dark and light distributions both have the 
same ellipticity. For practical purposes we fix the mass 
at M200 = 10 14 /i _1 Mq in this calculation, which corre- 
sponds approximately to clusters of the mean richness we 
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used (|Johnston et al.ll2007l . Table 6). The value used af- 
fects the concentration parameter and therefore the mass 
profile of the cluster, which affects the weighting of the 
map. If we use a value a factor of ten higher our light 
cllipticity results change by less than one sigma and in 
any case our main results, on the dark matter ellipticity, 
are changed imperceptibly because the uncertainties on 
those are dominated by the much larger uncertainty on 
the dark matter quantities. 

We calculate probabilities in the resulting three dimen- 
sional space (/lm, K and no) by calculating a \ 2 between 
the predicted number (Eq. I16|) and the observed number 



X 



E 



(n pre d(r l ) - riobs(ri)) 2 



(17) 



and calculating probabilities from this; Pr — e x / 2 . 

Our assumption is that the errors are Poisson, there- 
fore in the limit of large numbers the error on the 
number of galaxies, as used in the \ 2 calculation, is 
a rii = \J nP red (rj).We calculate our main results for each 
of the three stacked light maps shown in Fig. [TJ 

3.1.2. Estimation of the dark matter axis ratio (/dm) from 
the stacked shear map 

To estimate the ellipticity of the dark matter distribu- 
tion we calculate probability as a function of cluster axis 
ratio /dm and cluster mass. We calculate the probability 
from the \ 2 between the observed and predicted shears, 
using the uncertainty on the shear values from Eq. 1121 
We marginalize over the cluster mass with a flat prior to 
obtain the probability as a function of the 2D axis ratio 
/dm- 

3.2. Redshift distributions 

As discussed in Sec. [2j we divide our cluster sample 
up into four redshift bins. Due to the large photometric 
redshift uncertainties we decided not to use the redshift 
information in the shear galaxy catalogue. Therefore the 
'shear galaxies' may be in front of, behind, or part of 
the cluster. The shear for each lens-source pair depends 
on the redshift of both the lens (cluster) and the source 
(shear galaxy). To calculate our theoretical model, we 
need a prediction for the distance ratio in Eq. [14] at each 
possible redshift. This must be averaged, weighted by 
the number of galaxies at each redshift. In other words, 
we need to calculate 



\ D s 



f™(D ds /D s ) n(z s ) dz s 
Jo°° n ( z s) dz s 



(18) 



Note that the integration in the nominator starts at the 
lens redshift, zj,, so that galaxies between us and the lens 
do not contribute to the shear signal, as they are not 
influenced by the presence of the cluster. We est imate 
n(z s ) from Fig. 3 and Eq. 8 in lSheldon et all (|2001h . and 
as a result we use z c = 0.22. However our results are 
quite insensitive to these numbers because we focus only 
on the ellipticity of the dark matter halo and not on its 
mass. 

We obtain (D ds /D s ) = [0.51,0.37,0.26,0.17] for the 
four redshift bins. The values are low for high redshift 
bins because a large fraction of the galaxies are between 
us and the cluster, and therefore do not contribute to 



the shearing effect. The calculation of D ds /D s is ap- 
proximate because we assume all clusters to be located 
at the center of their redshift bin. However, this does 
not affect the ellipticity of the theory prediction. The 
distance ratio is incorporated into the predictions using 
Eq.[H 

3.3. Cluster decontamination 

Because spectroscopic redshifts are not available for all 
shear catalogue galaxies, there will always be a degree of 
contamination by cluster members in the shear signal 
(as seen in Fig. QJa)). Since we assume that the cluster 
members have no systematic alignment (but see Sec. [5] 
for an assessment of the implications of this assumption) , 
members that have leaked into the shear galaxy catalogue 
will tend to dilute the shear signal. We correct for this 
dilution in our analysis. The corrected shear is given by 

7cor = —7 19 
n 

where n(r) is the total number of galaxies in the shear 
galaxy catalogue a two dimensional position r from the 
center, no is the number of galaxies not in the cluster (see 
Eq. [T6|) and 7 is the observed, uncorrected shear. We use 
the best fit no value from the % 2 fit to the light matter 
distribution. Inside the cluster, we have n(r) > no, so 
the observed shear will be boosted by correcting for the 
contamination. 

We mask the central regions from our analysis for sev- 
eral reasons, (i) In the very central regions, cluster mem- 
bers will obscure the shear galaxies, (ii) The cluster cen- 
ter may be incorrect and thus the central parts may ap- 
pear erroneously circular after the stacking, (iii) There 
is actually an uncertainty associated with the correction 
factor that we have not taken into account. This error is 
due to Poisson fluctuations in the true number of non- 
cluster members per pixel, and will be more significant 
when the value of the correction factor is large, which oc- 
curs in the central region where the observed number of 
galaxies in the shear galaxy catalogue peaks sharply. For 
all these reasons, we mask out the central region using 
a circular mask with r mas k = O-5/i" 1 Mpc. The radius 
of the mask was set where the correction factor increases 
above 1.5 (as calculated for 0.20 < z < 0.25), shown 
by the dashed circle in Fig. [T] This corresponds to 3.3 
arcminutes at z = 0.225. 



4. RESULTS 

We now present our results on the stacked cluster el- 
lipticity, focussing first on the NFW profile. 

4.1. Results for the light matter axis ratio (/lm) 

The left panel of Fig. [3] shows the one-dimensional rel- 
ative probability of the axis ratio /lm, marginalised over 
K and no (see Eq. [T5]) . The solid lines represent re- 
sults from using the total cluster members, i.e. galax- 
ies in both the shear galaxy catalogue and cluster mem- 
ber catalogue (corresponding to the right hand panel of 
Fig [T]). The blue, dashed line shows results from us- 
ing only the (members-removed) shear catalogue galax- 
ies, corresponding to the left hand panel of Fig. [I] for the 
NFW model. This shows that the light matter distribu- 
tion is clearly elliptical with an axis ratio of /lm ~ 0.6. 



TABLE 1 

Axis ratio results for the light maps 



Population 


z 


SIE (ft m") 


NFW (ftiul 


Shear galaxies 
Total 
Members 


0.10 < 2 < 0.15 
0.10 < z < 0.15 
0.10 < z < 0.15 


0.548 + 0.013 - 0.013 
0.555 + 0.014 - 0.014 
0.302 + 0.006 - 0.006 


0.620 + 0.015 - 0.015 
0.614 + 0.012 - 0.013 
0.544 + 0.006 - 0.007 


Shear galaxies 
Total 
Members 


0.15 < z < 0.20 
0.15 < z < 0.20 
0.15 < z < 0.20 


0.672 + 0.011 - 0.011 
0.657 + 0.012 - 0.011 
0.261 + 0.005 - 0.005 


0.649 + 0.012 - 0.012 
0.631 + 0.010 - 0.010 
0.475 + 0.006 - 0.006 


Shear galaxies 
Total 
Members 


0.20 < z < 0.25 
0.20 < z < 0.25 
0.20 < z < 0.25 


0.581 + 0.012 - 0.012 
0.578 + 0.011 - 0.010 
0.254 + 0.004 - 0.004 


0.583 + 0.011 - 0.010 
0.575 + 0.009 - 0.009 
0.456 + 0.005 - 0.005 


Shear galaxies 
Total 
Members 


0.25 < z < 0.30 
0.25 < z < 0.30 
0.25 < 2 < 0.30 


0.637 + 0.010 - 0.009 
0.624 + 0.008 - 0.008 
0.235 + 0.003 - 0.003 


0.619 + 0.008 - 0.008 
0.602 + 0.007 - 0.007 
0.434 + 0.004 - 0.004 


Joint (total) 


0.10 < 2 < 0.30 


0.610 + 0.005 - 0.005 


0.602 + 0.004 - 0.005 



TABLE 2 

Axis ratio results for the dark matter distribution. 



2 


SIE (/ DM ) 


NFW (/dm) 


0.10 < 2 < 0.15 
0.15 < 2 < 0.20 
0.20 < 2 < 0.25 
0.25 < 2 < 0.30 


0.754 + 0.230 - 0.186 
0.522 + 0.159 - 0.115 
0.958 + 0.251 - 0.199 
0.853 + 0.306 - 0.215 


0.614 + 0.400 - 0.176 
0.269 + 0.158 - 0.054 
0.599 + 0.379 - 0.208 
0.614 + 0.444 - 0.266 


Joint analysis 


0.747 + 0.102 - 0.094 


0.480 + 0.136 - 0.086 



The errors on /lm are calculated by finding the 68% 
iso-probability limits from the probability distribution 
P(Ilm) after marginalizing over the other fit parame- 
ters. We find an error of ~ 0.005. 

Table Q] shows three light matter ellipticity results for 
each redshift bin and for each theoretical model (SIE and 
NFW): (i) for the shear catalogue galaxies with cluster 
members removed, (ii) for the total cluster members (iii) 
for the cluster members and (members from the cluster 
catalogue plus the extra members included in the shear 
galaxy catalogue). The results from the NFW and SIE 
are qualitatively similar. The last line in Table [T] shows 
the joint results combining all bins, for the total cluster 
members. 

There is a clear detection of ellipticity based on the 
number density of shear catalogue galaxies alone. We 
see that the distribution of galaxies in the members cat- 
alogue is more elliptical than that of the shear catalogue. 
This is not surprising because we stacked the galaxies ac- 
cording to the members catalogue. Even if the galaxies in 
the members catalogue had been sampled from a circu- 
lar distribution then the finite number of galaxies would 
provide a rotation direction so we would have in effect 
stacked the random noise to produce an ellipticity for the 
cluster members, while making the light map from the 
shear catalogue more circular. 

This extra induced ellipticity could be simulated, but 
in fact there is no need because the axis ratio from the 
shear galaxy catalogue alone is very similar to that based 
on the total catalogue containing both the shear cata- 
logue and the cluster member catalogue. Therefore we 
do not consider the results from only the cluster member 
catalogue any further. 

The similarity of the cllipticities from the shear galaxy 
catalogue and the total catalogue is largely due to the 



fact that there are very few galaxies from the members 
catalogue outside our mask radius. The important con- 
clusion for our paper is that any misalignment in the 
rotating and stacking will affect the light and dark mat- 
ter ellipticity the same. The ellipticity observed in the 
light map makes us believe that the rotation angles are 
not completely random and therefore that we can hope 
to detect some ellipticity in the dark matter, if indeed 
the dark matter halo is elliptical. 

4.2. Results for the dark matter axis ratio (7dm j 

Ellipticity results for the dark matter analysis can be 
found in Table [2j for each individual redshift bin as well 
as for the joint result combining all redshift bins. For the 
NFW model, the joint result of /dm = 0.48Io;Qg excludes 
a circular mass distribution (/dm = 1) by over 3er (if the 
probability distribution is Gaussian). The lowest axis 
ratio is in the second redshift bin (0.15 < z < 0.20). 
We re-analyzed redshift bin 2 by dividing it into sub- 
bins, but found the result unchanged on removing the 
clusters with the lowest axis ratio. The joint result for 
the remaining redshift bins was then /dm = 0.607^q 
This is a weaker detection than our final result including 
the second redshift bin, but still a tentative detection of 
dark matter ellipticity. 

The right panel of Fig. [3] shows the probability as a 
function of axis ratio for the dark matter distribution. 
The dark (red) line shows the result from using an NFW 
model, and the light (green) line shows the result from an 
SIE model. The vertical dotted line is showing /dm = 1 , 
which represents a circular mass distribution. 

We can see from the figure that the probability is not 
totally Gaussian, with a tail to larger /dm values. (Note 
that an axis ratio greater than unity corresponds to an 
ellipse aligned along the y axis.) We find that 99.6 per 
cent of the probability is below /dm = 1; therefore we 
consider this to be a reasonable detection of dark matter 
ellipticity. Note that the best fit NFW result is quite a bit 
more elliptical than the SIE result, but the NFW result 
has a longer tail extending to higher /oM-values than the 
SIE. For the SIE model 98.2 per cent of the probability 
is below /dm = 1- The NFW is a more realistic profile, 
so we trust these results most. 

Fig. H] shows 68% and 95% contours of the two- 
dimensional probability distribution for the axis ratio of 
the dark matter (/dm) versus the axis ratio of the light 
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Fig. 4. — 68% and 95% contours of the two-dimensional proba- 
bility distribution for the axis ratio of the dark matter (/DM)versus 
the axis ratio of the total galaxy number density distribution (/lmi 
representing the light matter). Light (green) contours show the re- 
sult from using an SIE profile as modelling the cluster, and dark 
(red) contours show result from using an NFW profile. The hor- 
izontal dotted line shows /dm = 1, and the dotted line rising to- 
wards the right shows /dm = /lm- 




DM LM 



Fig. 5. — One-dimensional probability distributions as a function 
of the ratio /dm//lm- Light (green) line is for the SIE result while 
dark (red) line is for the NFW result. 

matter (/lm)- The dotted lines also shown on the plot 
are /dm = 1 (horizontal dotted line) and /dm = /lm 
(dotted line rising towards the right). The two col- 
ors shown represent the different profiles used as cluster 
models: light (green) contours are for the SIE model, 
and dark (red) contours are for the NFW. We see that 
/dm = 1, that is, a circular dark matter distribution, 
crosses just within the outermost (95%) contour for the 
SIE model and the NFW contours are both below the 
/dm = 1 line. 

Fig. shows the one-dimensional probability distribu- 
tion as a function of the ratio of axis ratios of the dark 
and light matter: /dm//lm- Again, the light (green) 
line is for the SIE model profile and the dark (red) line 



is for the NFW model. The vertical dotted line shows 
where we have /dm = /lm, that is; the axis ratio as 
deduced from the shear values equals the axis ratio as 
deduced from the distribution of the cluster galaxies. 
Both profiles are consistent with having the same ellip- 
ticity in the light and in the dark matter, /dm//lm = 1. 
We find /dm//lm = 0.91±%H for the NFW model and 
/dm//lm = 1.26±g 1? for the SIE. The higher SIE result 
is driven by the larger /dm for the SIE, which makes the 
dark matter distribution appear more round. 

4.3. Dependence on close neighbor distance 

We made a decision to exclude all clusters with a more 
massive neighbor within an angular size corresponding to 
bh^ 1 Mpc at the cluster redshift. In this section we test 
how dependent our results are on this decision. Table [3] 
shows that our main result changes little as the cut is 
made 50 per cent smaller or larger. As expected, when 
a larger radius is used, fewer clusters survive in the cat- 
alogue, therefore the uncertainties become larger. Using 
a smaller radius does increase the ellipticity which is to 
be expected if the main effect is to include more physi- 
cally associated clusters which are more likely to lie along 
the major axis of the cluster, perhaps due to formation 
along an intervening filament. We find similar results for 
the light ellipticity. Our main results are unchanged by 
changing the cluster isolation criterion. 

4.4. Misalignment simulations 

In order to interpret our results we need to take into 
account possible misalignments during stacking. When 
we rotate clusters according to the distribution of the 
cluster members before stacking, we assume that the ori- 
entation of the light is correlated with the orientation of 
the mass. A result of /dm consistent with 1 could in- 
dicate either a circular mass distribution or a random 
alignment between the light used for stacking and the 
dark matter. In the latter case, the circular result would 
be caused by the stacking of many clusters with different 
misalignments between the dark and the light matter. 

In order to quantify this effect we simulate many clus- 
ters all with the same input axis ratio The mis- 
alignment angles between the light and the dark matter 
for the simulated clusters are random, with a standard 
deviation ag. For each value of / n and ag, we rotate 
and stack the clusters in the same way as for the SDSS 
data. Therefore if ag is large, the input ellipticity will 
be smeared out a great deal and output axis ratio / ou t 
will be close to unity. If ag is small the output ellipticity 
will be more similar (or in the case of zero misalignment, 
equal) to / in - 

Fig. [5] shows the output dark matter axis ratio as a 
function of degree of misalignment ag and input axis ratio 
/in- In order for the input ellipticity to be very elliptical 
(/ — » 0), the misalignment between the dark and light 
matter must be ~ 50° for our dark matter /dm of ~ 0.5 
(see Table [2J . We conclude that the misalignment angle 
must be less than ~ 50°. If there is any misalignment 
between the light and the dark matter, the dark matter 
will be even more elliptical than the results shown in the 
table and can be read off Fig. [6] for a given misalignment 
angle. 

5. CONCLUSIONS AND DISCUSSION 
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TABLE 3 

Effect of neighbor removal on measured ellipticity. Results are shown for the NFW profile only. 



Cut/ h- 1 Mpc 


No. clusters 


Light matter /lm 


Dark matter /dm 


2.5 

5 

7.5 


6934 
4281 
2542 


0.587 + 0.004 - 0.004 
0.602 + 0.004 - 0.005 
0.625 + 0.006 - 0.006 


0.459 + 0.155 - 0.077 
0.480 + 0.136 - 0.086 
0.614 + 0.205 - 0.139 
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Fig. 6. — Contours show output axis ratios for simulated mis- 
alignments between the light and dark matter. Axes show input 
axis ratio /j n and standard deviation ag (in degrees) of the angle 
distribution of the simulated clusters. 



We have used galaxy clusters from the catalogue 
of iKoester et all (l2007aD and shear maps as used in 
ISheldon et al.l (|2004[ ) to investigate the light and dark 
matter ellipticities of galaxy cluster halos. We rotate the 
selected clusters so that their major axes are aligned, and 
stack them according to the m ethod which is described in 
iNataraian fc Refregierl (|2000f) for galaxy-galaxy lensing. 
This is the first time that this method has been applied 
to cluster lensing. 

Through the pattern imprinted on the shear maps by 
intervening massive cluster halos, we have detected a 
dark matter ellipticity of these halos at a 99.6% level, 
with an axis ratio of /dm = 0.481°;^ from a joint \ 2 ~ 
analysis for an NFW model, using 4281 clusters between 
0.10 < z < 0.30. We have corrected for dilution of the 
shear signal by cluster members left in the shear galaxy 
catalogues but have masked out the central areas where 
this correction factor grows too large. 

The light matter distribution of the clusters, as traced 
by the number density of individual cluster members, 
is also clearly elliptical, with a joint axis ratio /lm = 



0.60 



+0.004 
-0.005 



for the NFW model. Using the shear cata- 
logue alone gives very similar results to using the shear 
catalogue concatenated with the cluster member cata- 
logue, which means that we are essentially comparing 
like with like when comparing dark and light map el- 
lipticities. This is because both the light matter (shear 
catalogue) and the dark matter have been stacked in the 
same way. 

We find that the ellipticity of the dark matter distribu- 
tion is consistent with the ellipticity of the galaxy number 
density distribution. Our result is limited by the uncer- 
tainty in the ellipticity of the dark matter distribution. 
The results for the NFW and SIE agree within the errors, 
but any differences could be attributed to a changing el- 



lipticity with radius. We have not attempted to measure 
this effect since the uncertainties are too large. 

Shear maps from neighboring clusters will influence 
each other, making the pattern more elliptical or less, 
depending on the redshift of the neighbors. In order 
to reduce this effect as much as possible, we use only 
the cluster with the highest number of members (N^®) 

when two (or more) clusters are closer together than 5/i _1 
Mpc, as seen on the sky. We find that increasing or de- 
creasing the minimum distance to a close neighbor by 50 
% does not significantly affect our main results. 

We have also simulated the effect that a possible mis- 
alignment between the light and dark matter could have, 
and concluded that, according to our results, the light 
and dark matter must be misaligned by less than ~ 50°. 

Our measurement is very insensitive to overall calibra- 
tion biases in shear measurements, since these are de- 
generate with cluster mass, rather than cluster elliptic- 
ity. It is unlikely that biases in shear measurement would 
vary with angle around the cluster: residual point spread 
function anisotropics would be oriented at random with 
respect to the cluster major axis and would cancel out 
on stacking; the cluster member light may leak into the 
postage stamps used to measure shears of background 
galaxies, but we remove the central region where the 
number of confirmed cluster members is significant. 

Possibly the biggest potential systematic is 
intrinsic alignment of clus ter members point- 
ing at the clust e r center (ICiqtti fc Duttal Il994t 



Kuhlen et alj l2007t iPereira et all 120081: iKnebe etal l 
20081 : iFaltenbacher et all l2008h which would be a 



problem for us because we have not made a significant 
attempt to remove cluster members from our analysis 
and our background catalogue is not particularly deep. 
A thorough assessment of this effect is beyond the 
scope of the current work. To first order this would 
make our observed shear maps less elliptical since the 
contamination by cluster members is greatest along the 
cluster major axis, and thus the strong gravitational 
shears expected along this axis will be partially cancelled 
out by the cluster members which have the opposite 
ellipticity since they tend to point at the cluster center. 
However a full assessment of this effect would also have 
to take into account the variation in intrinsic alignment 
with respect to the clu ster major axis, wh ich appears to 
be more complicated (jKuhlen et al.l l2007h and possibly 
weakens the degree of cancellation preferentially along 
the cluster major axis. 

In I Hopkins et al.1 (|2005h . the authors use a large-scale, 
high-resolution TV-body simulation to predict cluster el- 
lipticities and alignments in a ACDM universe. They 
find an ellipticity 



(e) = 1 - (/) = 0.33 + 0.05z. 



(20) 



for the redshift range < z < 3. This redshift evolution 
is negligible for the redshift range considered here, and 
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due to the large uncertainties we do not try to detect 
any trends with redshift. For a redshift in the middle of 
the redshift range of our cluster sample (z = 0.2), this 
formula yields an axis ratio of / = 0.66 which is in good 
agreement with our result of /dm = 0.48i ; 09 . 

IHq et all (|2006h have used numerical simulations of 
cluster formation with the aim of investigating the possi- 
bility of using cluster ellipticities as a cosmological probe. 
They find that the mean ellipticity of high mass clusters 
is approximated by 



e{z) = 0.245 



1 - 0.256 



0.9 



0.00246 



03 



(21) 



Using C8 (z) — 0.8 and Q m = 0.3, this gives an ellipticity 
of e ~ 0.2, or / ~ 0.8. This result is more circular than 
our main result. However, to make a proper comparison 
of our result with theory it would be necessary to make a 
theoretical prediction that takes into account our obser- 
vation method, especially the overlaps in the shear field 
due to close neighbors and the impact of selecting the 
most isolated clusters. 

iHoekstra et al.l (|2004f ) report on a first weak-lensing 
detection of the flattening of galaxy dark matter halos, 
using data from the Red-Sequence Cluster Survey. They 
find an average galaxy dark matter halo ellipticity of 
(e) = 0.33+q;q 9 . They also find a detection that dark 
matte r halos are rounder t han t he light map. In a later 
work, Mandclbaum et al. (2006) did not find a definite 
detec tion of this effect in the larger SDSS dataset. How- 
ever, iParker et all (|2007f ). measuring the ratio of tan- 
gential shear in regions close to the semi-minor versus 
that close to the semi-major axes of the lens, find some 
evidence of a halo ellipticity of ~ 0.3 using early data 
from the Canad a-France-Hawaii Tel e scope Legacy Sur- 
vey (CFHTLS). iMandelbaum et al.l (|2006h commented 
that the stronger signal in clusters of galaxies means 
that there is more chance of making a detection of el- 
lipticity in this higher mass data set. Therefore we have 
confirmed the detection of dark matter halo ellipticity, 
extending the measurement to cluster scales. We find no 
evidence for different ellipticities for the light and dark 
matter distribution o n clust er scales. 

IMandelbaum et all (|2006l ) focus on measuring the ra- 
tio of (dark matter) halo ellipticity to galaxy (light) el- 
lipticity £h/ e g where e = (a 2 — b 2 )/(a 2 + b 2 ). They find 



= 0.60 ± 0.38 (68 per cent confidence) for ellipti- 
cal galaxies, which is most comparable with our result 
of /dm//lm = 0-9lioi;5 (68 per cent confidence) where 
/ = b/a. Repeating our calculation using the different 
ellipticity parameters we find e^/e 9 = 1.37l ; 26 . It is 
expected that our value is greater than unity, since we 
find that the dark matter is more elliptical (lower axis 
ratio) than the light matt er. However it is no t signifi- 
cantly larger. The result of [Mandclb aum et al.l (|2006f ) is 
the opposite side of unity, but the difference is not sig- 
nificant ijJid_jvBjn^te_Jhat ourresult is for clusters and 
that of IMandelbaum et alj (|2006l ) is for galaxies. 

The anisotropy of the lensing signal around individ- 
ual galaxies, if definitely detected, has been said to 
pos e a serious problem f or alternative theories of grav- 
ity (|Hoekstra et al.ll2004l ). For galaxy clusters it is more 
complicated, as the dominant source of baryons is the in- 
tracluster gas. Lensing by clusters has been found to pose 
a problem for Modified Newtonian Dynamics (MOND), 
as there see ms to be a need to include a dark ma tter 
component (|Sandersll200llTakahashi fc Chiball2007f >. 

Because the gas - the dominant baryonic component 
of clusters - is collisional, we can suppose that on cluster 
scales it will be less elliptical than the light. We may 
end up concluding that either (i) the gas distribution is 
elliptical (which we would not expect) or (ii) our result 
is inconsistent with MOND. However, further study and 
simulations of this is clearly needed. 
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